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A procedure has been developed for computing flutter characteristics of single finite panels, 
particularly at low supersonic Mach numbers where static and quasistatic aerodynamic ap- 
proximations are not valid. Air forces from exact linearized potential-flow theory are used. 
The panel is considered as finely divided into many boxes, and the aerodynamic influence 
coefficients between all pairs of boxes are computed by numerical integration. The flutter 
analysis is a modal type of analysis which can be used with the box method for calculation of 
the flutter stability of any flat or nearly flat panel, whether of isotropic or anisotropic stiffness, 
and even for buckled panels for which the flutter is a small-amplitude, simple harmonic, 
superimposed motion to which linear theory would apply. Certain results are presented for 
flat unstressed isotropic panels. Variables whose effects were studied are panel lengtli-to- 
width ratio, Mach number, and the air cavity behind the panel. Panels with simply sup- 
ported and with clamped edges were studied. 


Nomenclature 

a. = speed of sound 

A v = aerodynamic influence coefficient giving the ve- 

locity potential at a box due to unit downwash 
on another box [Eqs. (3) and (9)1 
B, = number of boxes in stream direction 

B x , = number of boxes in cross-stream direction 

d = depth of cavity behind the panel 

g = structural damping coefficient of the panel 

hi,hj — distribution over the panel of mode shape for 

modes i and / [Eq. (14)] 

H,(x,y,r) = distribution of time-varying vertical deflection of 

panel for mode j [Eq. (14)] 
i = ( — 1 ) m = unit of imaginaries 

= integrals in A v [Eqs. (10-12)] 
k t = ue/V = reduced frequency with reference length e 
ki = ul/V = reduced frequency with reference length l 
l = length of panel in stream direction 

m,A = mass per unit of surface area of panel 

M = V /a, = Mach number 
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generalized mass for mode i 
Mi/lwmA (applicable for uniform panel only) 
number of modes in a flutter analysis 
perturbation pressure due to mode,/ [Eq. (17)] 
dynamic pressure of airstream 
time-varying generalized coordinate of motion for 
modes i and/ 

complex amplitude of q, [Eq. (14)] 
generalized aerodynamic force from the pressure 
due to mode / and the modal deflection of mode 
i [Eq. (18)] 

nondimensional computational quantity con- 
tained in Qa [Eq. (21)] 
thickness of panel 

transformed panel coordinates in x and y direc- 
tions, respectively, based on e/2 as a reference 
length [Eq. (5)] 

lower and upper limits of integration with respect 
to u and to v [Eq. (6)] 

downwash velocity at the panel surface for mode / 
width of panel in the cross-stream direction 
panel coordinates in stream and cross-stream 
directions, respectively (see Fig. 1), based on 
£ as a reference length in Eqs. (3) and (5) 
values of x and y at center of box m,n 
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Fig. 1 Plan view of panel divided into boxes with coordi- 
nate system, dimensions, and a forward-facing Mach cone 
with apex at box center # m ,y n * 
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| Subscripts 


* B 

*• 

s hi 


3 


m,n 


te 


panel coordinates based on reference lengths l 
and w, respectively 
see Eqa. (26-28) 

(A/s _ 1)1/2 
width of box 

dummy variable of integration for y 
panel-to-air mass density ratio 
dummy variable of integration for x 
density of airstream 
lB x ,/wB„ length-to-width ratio of box 
time 

velocity potential at center of box m,n 
velocity potential at center of box m,n due to unit 
down wash over box p,v 
frequency of flutter motion 

natural frequencies of panel modes i and j and of 
a chosen base or reference mode, respectively 
(w B /o>y( 1 + iff), flutter eigenvalue parameter 
M 2 k e /p 


base or reference value 
indices denoting mode numbers 
indices denoting streamwise and cross-stream 
numbered location of a box 
indices denoting relative locations of influencing 
and influenced box 
value at trailing edge of panel 
indices with same function as m,n 


Introduction 

T I HI ERE has been considerable interest and activity in the 
analysis of panel flutter. Perhaps the most recent re- 
views of the work done are those by Stocker 1 and by Fung 2 
which cover both theoretical and experimental investiga- 
tions. The present paper is concerned with an important 
type of panel and speed regime for which analysis has been 
lacking. That type of panel is the single rectangular one 
(embedded in a rigid surface), and the speed regime is the 
low supersonic, especially below Mach number 2 112 . For these 
Mach numbers, simple approximations for the aerodynamic 
forces, such as from Ackeret theory or piston theory, are 
not satisfactory. Use of the exact theory for linearized 
supersonic flow is a logical step. This exact theory was 
applied by Luke and St. John 3 to an effectively infinite- 
width panel separated into bays by supports equally spaced 
along the width. Zeijdel 4 obtained additional results for 


- ■ * 

such a panel with added supports that divide the length of 
the panel into equal bays. For this type of panel the general- 
ized aerodynamic forces needed in the flutter determinant 
can be obtained by evaluation of closed-form expressions. 
For finite-width panels, however, closed-form expressions 
have not been found for the aerodynamic force terms even 
for simple mode shapes of the panel. 

The subject of the present paper is the determination of 
the aerodynamic forces by, in effect, a numerical integration 
process and the subsequent use of the air forces to compute 
flutter boundaries. For this problem, current high-speed 
computing machines make practicable the obtaining of solu- 
tions by such “brute force” methods that would not be 
attempted otherwise. The panel is considered to be finely 
divided into a large number of equal-size boxes. The aero- 
dynamic influence coefficient of each box on each other box 
is calculated and used in a modal-type flutter stability 
analysis. The mode shape properties can be supplied from 
either calculation or experimental measurement. 

A number of flutter results are presented for flat unstressed 
panels with isotropic stiffness and by use of calculated mode 
shapes. Calculated modal frequencies were used, with one 
exception. The principal variables studied are the panel 
length-to-width ratio and the flow Mach number. Effects 
of the air cavity behind the panels also are examined in a 
qualitative sense. 

Although the results presented are limited to simple panels, 
this modal analysis with a box method has a broad applica- 
tion and lends itself to the calculation of the flutter stability 
of any flat or nearly flat panel, whether unstressed or stressed 
(as by thermal expansion) , whether of isotropic or anisotropic ' 
stiffness, and even for buckled panels for which the flutter 
is a small-amplitude, simple harmonic, superimposed motion * 
to which linear theory would apply. 


Analysis 

The panel to be analyzed and the coordinate system are 
shown in Fig. 1. It is a single rectangular panel of length 
l, width w. and the side edges are aligned with the remote 
wind direction. It is a flexible part of an otherwise inflexible 
surface which extends at least far enough so that upper and 
lower surfaces of the panel induce no aerodynamic effects 
upon each other. The supersonic stream passes over the 
upper surface. No account is taken of any perturbation 
pressures on the panel lower surface except for some re- 
marks regarding cavity effects near the end of the paper. 

The panel is divided into a number of equal-size rectangu- 
lar boxes. The number of boxes in the stream and cross- 
stream direction are B. and B x „ respectively. For con- 
venience and for purposes of reference in the computational 
procedure, the boxes are numbered in sequence, beginning at 
the box nearest the origin of coordinates. In the stream 
direction the index m = 0, 1, 2, . . . B, — 1, and in the cross- 
stream direction n = 0, 1, 2, . . . B x , — 1. A second set of 
box index numbers is needed, and these are \j. = 0, 1 , 2, . . . 
5,-1 and v - 0, 1, 2, . . . B x , — 1. Thus, the aerodynamic 
influence upon any box m,n due to any other box p,v can be 
referred to. 


Aerodynamic Influence Coefficients 

It is assumed that the panel is divided into a sufficiently 
fine gridwork of boxes so that the downwash over any one 
box can be taken as uniformly distributed at any instant 
and that the resulting pressure perturbation at the center 
of each box is a sufficiently accurate average of the pressure 
distribution over that box. 

For convenience in the computational procedure, a refer- 
ence length is chosen as e, the width of a box; that is, e = 
w/B x „ The velocity potential at the center of box m,n due 
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to a uniformly distributed downwash w(y,v) over the box 
y,v can be expressed for a simple-harmonic time variation as 

<p(m,n) = ew(n,v)A v (r,s) - ' (1) 

where 

r = m — y s = n — v (2) 

and, from Eq. (34) of Ref. 5, 


both sides of the Mach cone). For this condition 

7 " - \L ^ m/2)uj {^) du (io) 

The second form of Eq. (9) occurs for portions of a box cut 
by one side of the Mach cone (v x n), so that the limit i> 2 = 
'u, Vj = 2s — 1 & 1. Since t> 2 = u. then cos -1 (t> 2 /u) = 0 and 


.4.M . Iff 


C 0 s(fi/M) [(x m - 0 2 - PKVn - y) 2 ] 112 
„ Re[(x m -~£) 2 - 0*(y„ - y) 2 ] 112 


d£ di i 


(3) 


in which the area of integration S^„ is the area of box y,v. 
In the denominator, the indicated real part of the radical 
term has a nonzero value for any part of .S Mi „ which lies ahead 
of the forward facing Mach cone with apex at x„, y„ (see 
Fig. 1). All boxes and parts of boxes lying behind the Mach 
cone have a zero contribution to 4,(r,s). The frequency 
Mach number parameter 0 is 

12 = M% c /P (4) 

where the reduced frequency k e = coe/F contains the chosen 
reference length e. 

As in Ref. 5, a transformation of variable is made as 
follows : 


Im 


« = e cos-i ill 



cos(j3fl/27l/)(M 2 - t> 2 ) 112 - 1 
yu- — v-)'- r - 


dv)du 

/ 


( 11 ) 


The third form I G s occurs for boxes that are completely ahead 
of the Mach cone and also for portions of boxes ahead of the 
point where the Mach cone passes out through the side of the 
box: 


Ig 3 = ~ f e W*>*( cos- 1 — — cos -1 — + 
2n*' u » \ u u 


Then 


x m — £ = fiu/2 y„ — jj = v/2 


(5) 


T-cos((3S2/2M)(m 2 - v *) 112 - 1 , \ , 

J VI ( w 2 _ v t)m dv ) du 


(12) 





e-i(pa/ 2 )u du X 



cos(;8ft/2.U)(w 2 - v 2 ) 1 ' 2 
(m 2 - v 2 ) w dV 


(6) 


where the upper and lower limits of the surface integrals, 
indicated by v u t> 2 , u h and w 2 , are determined for any given 
box by its edges except where the forward Mach cone emanat- 
ing from x m , y„ passes through the box. 

The integral of Eq. (6) is the one to be evaluated nu- 
merically. It is to be noted that the integrand contains a 
singularity of order (— ^) at u = v, that is, along the Mach 
lines from x„, In order to avoid the potential inaccuracies 
involved in integrating in the vicinity of this singularity, 
the numerator of the integrand is replaced as follows: 


cos(/?f2/2./lf)(w 2 — f 2 ) 1/2 = 

[cos(0O/2i¥)(w 2 - v 2 ) 112 - 1] + 1 (7) 


The integration with respect to v can be partly carried out 
formally : 



cos03R/2.1/)(w 2 - v 2 ) 1 ' 2 - 1 
(w 2 - a 2 ) 1 ' 2 


dv 


i; r m cos(|8SV24/)(w 2 — 

„ ' I 77 . 2 _ .,2V 


„2U/2 


Table 1 shows all possible relative locations of a box and the 
Mach cone. The applicable integral, whether 7 0l , / 02 , or 
1 (33 is shown, and the limits of integration v u i> 2 , u u and u, 
are listed. The limits shown are consistent with the refer- 
ence length e chosen for Eq. (1) and with the transformation of 
Eq. (5). The quantity <r = lB„/wB, is the length-to-width 
ratio of a box. 

With the assumptions that have been made, the aero- 
dynamic influence coefficients have right-left symmetry 
about s = 0; that is, 

A^(r,s) = A v (r,-s) (13) 

Thus, all of the A v (r,s) which are needed for any box are 
obtained for either one of the rear corner boxes. 

Once all A„(r,s) are computed, the total <p,(m,n) for any 
downwash mode j at the center of any single box m,n is ob- 
tained by a matrix multiplication : 

<Pi(rn,n) = Fe^E^^) 

H v ' 

= Ve |^^J{/l v (r, S )) (14) 

in which any element A v is zero where m — y = r would be 
negative. 

Downwash 


The integrand of the remaining integral with respect to v is 
zero at u = v, and the integral can be evaluated more ac- 
curately. Substitution of Eq. (8) into A v (r,s) gives the re- 
sult 


-4„(r,s) 


,s) = ~ f e-'^/^fcos- 
2? tJ u\ \ 


— cos 


-iii? 


+ 


f' 

U Vi 


! cos(|80/2M)(m 2 — v 2 ) 112 — 1 


( m 2 - v 2 ) 


1/2 


dv ) du (9) 


There is a need for only three forms of Eq. (9). The first is 
for the condition that v t = —v? = — u (which can occur only 
for v = n, and then only for the portion of a box y, v cut by- 


The downwash mi, is employed as follows. Let the panel 
deflection in mode j be expressed for simple-harmonic time 
variation as 

Hj(x,y,T) = qMhfay) = qi(T)e iwT hj(x,y) (15) 

where 7, is the complex amplitude of the generalized coordi- 
nate of motion, and hj(x,y ) is the shape of the deflection mode. 
The downwash ratio is 


Wj 

V 


£>//, 1 dH, 

b(lx) + V br 




(16) 
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where x is a convenient streamwise coordinate that is zero at 
the panel leading edge and 1.0 at the trailing edge. Substi- 
tution of Eq. (16) into Eq. (14) gives the result 

= Vt j + i ^A, j{A„(r,s)} (17) 

in which A, is hj(x„y v ). Thus, the <pj(m,n) are obtainable in 
a systematic manner for each mode j used in a flutter analysis. 


flutter determinant is the integral over the panel surface of 
the product of the deflection mode shape hi times Ap,-; that is, 

Qa = ~ jff hi (^r + i y d(lx)d(wy) (19) 

The coordinate y ranges from zero at the left edge to 1.0 at 
the right edge. If the term hi(b<pj/bx) is integrated by parts, 
there is obtained 


Generalized Aerodynamic Forces 

The perturbation pressure Ap, is obtained from the 
velocity potential by 


A Pi = —pi 


y <><pj b<Pf ~ 

L o(ix) or _ 


- + >#*>,) OS) 

i \ Ox v / 


Qa = - y f^hi(x u ,y)ip,(x u ,y) - 

f <pi ~ i v ] < ^ w£) ( 2 °) 

Since hi(x,„y) is zero for panels supported at the trailing edge, 
On = f f (Pi ( - t - h.\ d(li)d(wv) (21) 

l JJ ■■ \ bx V / ' ' ' 


The generalized aerodynamic force term Qi, needed in the 

Table 1 Types of integrals and limits of integration for 
computing A v (r,s) for all possible relative locations of the 
box y, v and the Mach cone from box m,n 


Limits of integration 
V 1 v 2 U 1 “2 




Limits of integration 


2s-l u (2r-l)o/p (2r+l)o/B 


In Eq. (21) the quantity [(M./dx) — i(ul/V)hi] readily is 
seen to be the complex conjugate of the amplitude of the 
downwash ratio for mode i [where the subscript t is not to 
be confused with the unit of imaginaries i — (— 1) I,J ] . 

In practice the integral of Eq. (21) is evaluated by a sum- 
mation or matrix multiplication that can be put in the con- 
venient form 


Qii = pVi Bjh)U Qii * 


( 22 ) 


where 


Ur _< £*■]{[» +*‘r*']i' l 4 <23) 


Solution of the Flutter Determinant 

For the present purpose it will be assumed that the modes 
to be used in the flutter solution are orthogonal, that is, 
with neither inertia nor elastic coupling between modes. 
The equation expressing the dynamic equilibrium of motion, 
derived from application of the Galerkin method, is 

[co 2 - W> 2 (1 + ig)]Miqi + = 0 (24) 

j 

(t = 1,2, . . . N) 

where w; is the natural frequency of mode t, N is the number 
of modes used in the analysis, g is an assumed damping coeffi- 
cient used as a flutter indicator, and Mi is the generalized 
mass of mode i: 

Mi — Iw / 0 J 0 [hi{x,y))Hx dy (25) 

with m A representing the panel mass per unit area. 

With use of the definition 0 = (co B /a;) 2 (l + ig) as the 
flutter eigenvalue parameter, where co B is any chosen base 
or reference frequency, Eq. (24) can be rearranged in a form 
convenient for eigensolution : 



[ -« + (2r)1 * + (5) K§r. ^ - 0 » 

where 



and k t = c ol/V is the reduced frequency based on the reference 
length l. For panels with a uniform mass distribution m A , 
the quantity a/ Mi in Eq. (26) can be replaced as follows: 

a/ Mi = (a*/M*)(l/p) 


(27) 
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where 


a* = v)/lB,B x r 

n i ( 28 ) 

, hm d y 

and 1 /m = pl/m A is the air-to-panel mass ratio. Thus, one 
now has all the necessary quantities for computing flutter 
boundaries. Panel dimensions and natural frequencies are 
employed only as ratios. The downwash ratios and other 
parameters involving the mode shapes are in forms that are 
independent of the panel size. 

As is commonly the case with exact air forces, the flutter 
boundary' cannot be computed directly, and an indirect 
method of computing and cross-plotting is necessary. In 
advance of a flutter solution for any given panel and Mach 
number, choices are made for B„ B z „ and the number of 
modes in the analysis. The downwash quantities bhi/bx 
and hi are computed or obtained from experiment at the 
center of each box for each mode i (or j) and arranged sys- 
tematically for use in the matrix multiplications. The quanti- 
ties Mi and a (or Mi*, 1/m, and a* for panels with uniform 
Ota) are computed, except that allowance for varying the air 
density in a (or 1 /m) is retained. 

Following judicious choice of the reduced frequency ki, 
the block of values of A v (r,s), then the <Pi(m,n ) at each box 
center for each mode, and then the generalized aerodynamic 
forces Qa are computed and employed in the flutter deter- 
minant. For each of a number of values of the air density, 
the complex eigenvalues 0 are computed. There is one 
eigenvalue corresponding to each chosen panel mode. By 
plotting curves of g against 1 /m for a sufficient range of air 
densities, the existence or nonexistence of an eigenvalue with 
g = 0 for some air density is established for each mode. An 
associated curve of the stiffness parameter l/V or oj 1 l/a, 
is plotted against mass ratio — one curve for each eigenvalue. 
(The use of the speed of sound a, rather than V permits direct 
comparison of flutter boundaries for various Mach numbers.) 

On a plot of 1 /m against stiffness parameter, points now 
can be placed which represent g = 0 or some other small value 
of g believed representative of the actual hysteresis structural 
damping of the panel. Each point is on a separate stability 
boundary. 

By successive judicious choices of hi, a series of points can 
be placed on the plot of 1 /m against stiffness parameter and 
stability boundaries drawn through the points. Four differ- 
ent types of stability boundaries have been found, in general, 
and three of them are illustrated in Fig. 2. For each type, 
a boundary for g = 0 and one for g = 0.01 (as illustrative of 
some small value) are shown, and the arrows point in the 
direction of increasing reduced frequency ki. The type of 
boundary of Fig. 2a predicts that, if g were zero, a fixed panel 
thickness would be needed as the air density tends toward 
zero, but if g is not zero, panel thickness goes toward zero as 
air density does. 

For the type of boundary of Fig. 2b, the unstable region 
shrinks with increasing g and vanishes for some value that 
may be extremely small or as large as g — 0.05 or possibly 
larger, depending upon the panel and flow parameters. The 
crossing point of the boundary for g = 0 in both Figs. 2a 
and 2b through 1/m = 0 occurs at a value of ki for which the 
imaginary part of Qn(k t ) passes through zero. For the type 
of boundary of Fig. 2c, small increases of g have little or no 
effect, and this small effect can be either stabilizing or de- 
stabilizing, depending on the particular panel and flow param- 
eters. An important point regarding the boundary of Fig. 
2c, as well as the dashed boundary of Fig. 2a, concerns their 
resemblance to a parabola. For any boundary or portion of 
a boundary which is a parabola accurately described by 
1 /m = (wi l/a,) 1 times a constant, the air density and the 
panel thickness ratio t/l to prevent flutter are related by the 
formula </Z(p) 1/3 = const. Such a relationship is contained 





Fig. 2 Three general types of stability boundaries for g = 
0 and for g = 0.01. Arrows indicate direction of increasing 
reduced frequency k|. 

in the panel flutter parameter ( E(3/q) il3 (t/l ), which has been 
evolved and used by several investigators. The fourth 
possible type of boundary found from a flutter determinant 
is one that falls entirely in the negative I/m region, with 
the positive 1/m region being stable with respect to it. 

For any given panel and flow parameters, a set of stability 
boundaries is computed according to the number of modes 
used in the analysis. The flutter boundary separates the 
region that is stable with respect to all stability boundaries 
from the region that is unstable with respect to at least one 
stability boundary. It is the duty of the investigator to 
establish that the flutter boundary is converged, that is, that 
enough modes have been used in the analysis so that addi- 
tional modes do not alter the flutter boundary in any im- 
portant respect. 


Results and Discussion 
Effect of Panel Length-to-Width Ratio on Flutter 

Results have been calculated for several panel lengt.h-to- 
width ratios in the range from l/w = 0 (the two-dimensional 
panel) to 10, for M = 1.3, and with g = 0. The panel edges 
were either clamped or simply supported (pinned). The 
results for an aluminum panel at sea level are shown on Fig. 
3. The abscissa is the length-to- width ratio, and the ordinate 
is a panel flutter parameter ( fiE/q) 1,3 (t/l ). Analytical re- 
sults calculated for M — 1.3 and with g= 0 are shown for 
panels with clamped edges by the square points connected 
by the long-dash line and for panels with simply supported 
(pinned) edges by the round points connected by the short- 
dash line. Also shown for comparison on Fig. 3 is the ex- 
perimental flutter envelope from NASA TN D-451 6 which was 
obtained from tests on many panels with different sizes and 
materials, for Mach numbers from 1 .6 to 6, at different tem- 
peratures, and with various compressive loadings and amounts 
of buckling. The analytical results fall below the experi- 
mental envelope as would be expected, since the analysis 
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Fig. 3 Analytical variation of panel flutter parameter with 
1 / tv for aluminum panels at sea level compared with the ex- 
perimental envelope of Ref. 6. Present results for a single 
panel and previous results for an infinite span array of 
panels are shown for M = 1.3 and g = 0. 


is of flat unstressed panels, and the experimental envelope is 
established largely by marginally buckled panels. Also 
shown as a matter of interest are a few results from Refs. 3 
and 4 for an infinite span array of panel bays. 

A point of interest concerns the flutter frequency and 
flutter mode. For all the length-width ratios at least up 
through 4 on Fig. 3, the flutter frequency was near and usu- 
ally slightly above the first-natural-mode frequency. But 
even where the flutter frequency was almost exactly equal 
to the first natural frequency, the flutter mode was a coupling 
of more than one natural mode. The degree of coupling was 
small for the low values of l/w but increased with increasing 
l/w. For l/w = 10, there was a great change in the flutter 
mode. Twelve modes were required in order to obtain a 
converged flutter boundary. The flutter frequency was 
between the eighth and ninth natural frequencies. The 
largest component present in the flutter mode was the ninth 
natural mode, and the proportions of modes near to the ninth 
were also large. 


Effect of Mach Number 


be “two-dimensional.” In Fig. 21 of Ref. 7, the experi- 
mental thickness ratio to prevent flutter is shown as a func- 
tion of Mach number for a condition of no pressure differ- 
ence through the panel. The increase of thickness as M 
decreases below 2 I/2 was found to be much smaller than 
predicted by analysis of a two-dimensional panel. Even 
though the panel motion showed no cross-stream variation, 
it seems evident that, the aerodynamic forces exerted on the 
panel are not truly two-dimensional, especially near the panel 
side edges. The reason is that from the panel side edge to 
the tunnel side wall there is an area of nonmoving tunnel 
surface of width equal to about 6.2% of the panel width. If 
the tunnel side wall acted as a reflector, the image of the 
panel beyond the side wall would be at a distance of 12.4% of 
panel width from the panel. 

Since the images could not be accounted for with the 
present analysis, a single square panel was analyzed as a 
matter of interest. Six modes were used for the panel 
clamped at leading edge, simply supported at trailing edge, 
and free at side edges. The measured frequencies for the 
first two modes (the only ones available) and calculated fre- 
quencies for modes 3 to 6 were used. Flutter calculations 
were made at three Mach numbers: 1.16, 1.20, and 1.30. 
The results are shown on Fig. 4 for comparison with the ex- 
perimental results from Ref. 7. The upper curve is for g = 0, 
the next lower one for g = 0.01, and below that one for g = 
0.015. Thus, in terms of thickness ratio the analytical re- 
sults with a small value of g agree rather well with the ex- 
periments. However, there is a significant difference in the 
flutter frequency. The experimental frequency was below 
the first natural frequency and within about 15% of it. 
The analytical frequency for a single panel was slightly below 
the second natural frequency for M — 1.16 and 1.20 and was 
slightly below the third natural frequency for M = 1.30. 
Obviously, some factor or factors are not accounted for in 
the analysis to match the experiments, whether it is the 
boundary layer over the panel (as suggested in Ref. 7), the 
reflected images of the panel through the boundary layer on 
the tunnel side walls (see sketch in Fig. 4), the air in the 
cavity behind the panel, or the air puffing in and out of the 
cracks along the panel side edges. 


The thickness ratio t/l and the panel flutter parameter 
have been computed for an aluminum panel at sea level 
with l/w = 2, edges clamped, and for several Mach num- 
bers from 1.1 to 2.0. The results are given in Table 2. The 
thickness ratio to prevent flutter is virtually constant for 
M of 1.2 to 1.5, is 5% higher at M = 2.0, and is 6% lower at 
M = 1.1. This trend for low supersonic Mach numbers is in 
sharp contrast to the result for two-dimensional panels for 
which a great increase in thickness ratio is predicted to be 
required for Mach numbers less than about 2 1/2 . The ratios 
of flutter frequency to first natural frequency fell between 
1.0 and 1.1 for this panel and Mach number range. 

Reference 7 (Lock and Fung) reports the results of panel 
flutter experiments at low supersonic Mach numbers from 
1.16 to 1.45. The panels tested were square or nearly so 
(l/w = 0.95, 1.00, and 1.06), were clamped at their leading 
edges, approximately simply supported at their trailing edges, 
and free on their side edges to allow the flutter motion to 


Table 2 Clamped-edge aluminum panels at sea level 


M 

t/l 

(, &E/qY'\t/l ) 

1.1 

0.00292 

0.212 

1.2 

0.00311 

0.241 

1.3 

0.00312 

0.251 

1.4 

0.00311 

0.248 

1.5 

0.00311 

0.247 

2.0 

0.00327 

0.248 


Effect of Air Cavity behind the Panel 

Panels on aircraft and missiles generally have some kind of 
a space behind them which can have an effect on the flutter 
stability of the panel. Such an effect was reported in NASA 
TN D-827 8 on the flutter of a Fiberglas sandwich panel with 
a plastic-foam core. The panel had a length of 33.38 in. 
and a width of 20.31 in. Most of the tests were made with 
cavity depth of 1.5 in. behind the panel, and flutter was ob- 
tained over a Mach number range from 1.76 to 2.87. When 
the cavity was filled with layers of plywood so that the back 
surface of the panel touched the plywood, no flutter was ob- 
tained even though the airstream dynamic pressure was 
doubled. But when the cavity depth was made 0.5 in., 
flutter occurred at a dynamic pressure 40% lower than with 
the 1.5-in. cavity depth. 

Since the flutter frequency was much lower than the lowest 
acoustic-resonance frequency of the cavity, a simplified study 
was made on the basis that the only action of the air in the 
cavity is as an air spring that resists panel deflection in any 
mode that would compress the air but does not affect modes 
not involving compression. In the analysis the odd-num- 
bered modes (1, 3, 5, etc.) involve compression, and the even- 
numbered modes (2, 4, etc.) do not. The air-spring effect 
is largest on mode 1, much less on mode 3, and essentially 
negligible on higher modes. The analytically determined 
effect for a progressively decreased cavity depth primarily is 
to force the first-mode frequency upward toward the un- 
affected second-mode frequency. For sufficiently small 
cavity depth, the first-mode frequency could coincide with 
or even go above the second-mode frequency. One finding 
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Fig. 4 Comparison 
of analytical and ex- 
perimental t/l as 
functions of M for a 
square brass panel 
with clamped lead- 
ing edge and simply 
supported trailing 
edge. 


of the analysis was that determination of the natural fre- 
quencies of the panel in place over the cavity should be meas- 
ured with the air density equal to (or over a range of density 
including) that in the cavity during the flutter experiment. 
This is true because the test density is often a small fraction 
of the density of the outside atmosphere (where the still-air 
panel modes and frequencies usually are measured), and be- 
cause the incremental effect of the cavity air spring is propor- 
tional to the cavity air density. 

Since experimental frequencies were not available from 
Ref. 8 for the air densities of the flutter tests, flutter results 
have not been computed as a function of cavity depth as a 
variable. Instead, the panel flutter parameter has been 
obtained as a function of an assumed variable first-mode 
frequency. The result is shown qualitatively in Fig. 5 as 
the panel flutter parameter vs the first-mode frequency. 
As the first-mode frequency increases, the flutter parameter 
increases at an accelerating rate, reaches a maximum value 
in the vicinity of coi = W 2 , and then decreases sharply with 
further increase of o»i. It seems clear that the same trend 
would be obtained for decreasing cavity depth. 

Finally, there are undoubtedly other effects of the cavity 
which would become important for sufficiently small cavity 
depths. These are the virtual inertia and the frictional re- 
sistance of the air being pushed back and forth in the cavity 
by panel modes higher than the first. There could also be 
acoustic resonance effects at sufficiently high frequencies. 
The effect of the virtual inertia of the cavity air would be to 
lower natural-mode frequencies. In Ref. 8 it is reported 
that decreasing the cavity depth from 1.5 to 0.5 in. lowered 
the frequencies of modes 2, 3, and 4 by appreciable amounts, 
but the first-mode frequency could not be found. These fre- 
quency measurements were made at essentially sea-level air 
density for which the cavity effects are proportionately larger 
than for the much lower air density of the flutter test. 

Concluding Remarks 

A panel flutter analysis has been developed for which the 
generalized aerodynamic forces are computed by the “brute 
force” technique of dividing the panel into a larger number of 
boxes and computing and employing the aerodynamic in- 
fluence of each box on each other box. Such a technique 
can be resorted to where the aerodynamic forces are not 
obtainable from simple or elegant closed-form mathematical 
expressions. An important area of application of this 
technique is to finite panels at low supersonic Mach num- 
bers. Either experimental or analytical mode shapes and 
frequencies can be used. Thus, this type of analysis can be 
applied to any flat or nearly flat panel, whether unstressed 
or stressed (as by thermal expansion), whether of isotropic 
or anisotropic stiffness, and even to a small-amplitude flutter 
superimposed on a buckled deflection. 

Results are presented only for flat unstressed isotropic 
rectangular panels with side edges parallel to the air-stream 
direction and for a condition of no pressure difference that 
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Fig. 5 Qualitative effect on t/l and panel flutter param- 
eter of an increase of the panel first-mode frequency, 
such as could be caused by a decrease of the cavity depth 
behind the panel. 

would tend to bulge the panel. A panel flutter parameter, 
involving panel stiffness and stream dynamic pressure, is 
plotted as a function of panel length-to-width ratio for both 
clamped-edge and simply supported panels and is compared 
with the experimental flutter envelope from Ref. 6. In 
terms of panel thickness to prevent flutter, the analytical 
results fall moderately below the experimental envelope, 
which result is to be expected because the envelope was es- 
tablished largely by the more severe condition of marginal 
thermal buckling. 

The effect of variable Mach number was studied for a 
clamped-edge aluminum panel at sea level with a length-to- 
width ratio of 2. The thickness ratio was essentially con- 
stant for M = 1.2, 1.3, 1.4, and 1.5, was 5% higher at M = 
2.0, and was 6% lower at M = 1.1. Another variable Mach 
number study was made to compare with a series of experi- 
ments in a wind tunnel on square and nearly square panels 
at Mach numbers from 1.16 to 1.45. Analytical results 
agreed well with the experimental for the thickness ratio 
but not for the flutter frequency. 

The effect of a cavity behind the panel was analyzed quali- 
tatively on the basis of the cavity acting solely as an air 
spring to resist panel modes that would compress the air, 
thereby raising those modal frequencies. If the ratio of 
cavity air density to cavity depth is high enough, the first- 
mode frequency can be raised to equal or even exceed the 
second-mode frequency. An equality of first- and second- 
mode frequencies is generally a very unfavorable factor and 
tends strongly to produce flutter. 
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